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Abstract 

Within the Excited Baryon Analysis Center we have performed a dynamical coupled-channels analysis of the 
available p{e,e'Tv)N data in the region oi W < 1.6 GeV and < 1.45 (GeV/c)^. The channels included are 
7*A'^, ttN, rjN, and mvN which has ttA, pN, and aN components. With the hadronic parameters of the model 
determined in our previous investigations of nN ttN reaction, we have found that the available data in the 
considered W < 1.6 GeV region can be fitted well by only adjusting the bare N* helicity amplitudes 

for the lowest N* states in P33, Pn, Sn and D13 partial waves. The meson cloud effect, as required by the 
unitarity conditions, on the -y* N —> N* form factors are examined. 
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1 Introduction 



The electromagnetic parameters characterizing 
the excited nucleons, N* , and in particular the 
7* TV — > N* form factors, are important informa- 
tion for understanding the hadron structure within 
Quantum Chromodynamics (QCD). Thanks to the 
efforts in recent years, as reviewed in Ref. the 
world data of 'y*N — > A(1232) form factors are con- 
sidered along with the electromagnetic nucleon form 
factors as the benchmark data for developing hadron 
structure models and testing predictions from Lattice 
QCD calculations (LQCD). 

In the present work, we explore the extent to 
which the available p{e,e'Tr)N data in < 1.6 GeV 
can be used to extract the 7* — > iV* form factors for 
the N* states up to the so-called "second" resonance 
region. 

We employ a dynamical coupled-channels model 
developed in Refs. and extend our analy- 

sis of pion photoproduction reactions. We therefore 
will only recall equations which are relevant to the 
coupled-channels calculations of p(e, e'7r)A^ cross sec- 
tions. In the helicity-LSJ mixed-representation where 
the initial state is specified by its helicities 
and Ajv and the final MB states by the {LS)J an- 
gular momentum variables, the reaction amplitude of 



J* (q,Q^) + N{—q) ^ ^{k) + N{—k) at invariant mass 
W and momentum transfer = ^9'''?^ = — lj'^ 
can be written within a Hamiltonian formulation ^ 
as (suppress the isospin quantum numbers) 



LS N ,X-y X 



(k,q,W,Q')^ti 



R,J 



ik,q,W,Q^) 



+tLL.N.,^,Jk,q,W,Q') (1) 



where S'jv = 1/2 is the nucleon spin, W — uj+Ej^{q) is 
the invariant mass of the 7*A'^ system, and the non- 
resonant amplitude is 



{k,q,Q') 

+ ^ ^ / k''^d,k'tlg^^j^j^,^,j^j,g,{k,k',W) 

M'B' L'S' 

xGm'b' (fc', W)vi,s,M's'.x,x, {k',q, Q') (2) 

In the above equation, GniiB'{k' ,W) are the meson- 
baryon propagators for the channels M'B' = 
TTN,riN,TTA,pN,(TN. The matrix elements 

«/sMs,A-,A„('=>9>Q^). which describe the 7iV -> MB 
transitions, are given explicitly in Appendix F 
of Ref. The hadronic non-resonant ampli- 

tudes t-[sN^N,L's'M'B'ik,k' ,W) are generated from 
the model constructed from analyzing the data of 
TT nN reactions ^ . 
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The resonant amplitude in Eq. JT]) is 



D,AW)fi.,,^,Jq,W,Q''), (3) 

where the dressed N* ^nN vertex ^N^^LS^nNi^jW) 
and N* propagator Di,j{W) have been determined 
and given exphcitly in Ref. '4]. The quantity rele- 
vant to our later discussions is the dressed 7*A^ N* 
vertex function defined by 



photon 3-momentum (f^ in the laboratory frame as 
e=[H-2(|fe|VQ2)tan2(6le/2)]-i; K is the helicity of 
the incoming electron; (/>* is the angle between the 
TT-N plane and the plane of the incoming and out- 
going electrons. The quantities associated with the 
electrons are defined in the laboratory frame. On the 
other hand, structure functions of "f*N ttN pro- 
cess, = c^„(W^,Q^cos6';) (a = T,L,LT,TT,LT'), 
are defined in the final ttN center of mass system. 
The formula for calculating aa from the amplitudes 
defined by Eqs. (P)-© are given in Ref. 

In this first-stage investigation, we only consider 
the data of structure functions of p(e, e'7r°)p tiS'^ 



(4) 



Similar to what was defined in Ref. 0, S] j we call 
the contribution of the second term of Eq. ^ the 
meson cloud effect to define precisely what will be 
presented in this paper. We emphasize here that the 
meson cloud term in Eq. (jj]) is the necessary conse- 
quence of the unitarity conditions. How this term 
and the assumed bare N* states are interpreted is 
obviously model dependent. 

Within the one-photon exchange approximation, 
the differential cross sections of pion electroproduc- 
tion can be written as 
da'' 



dE,,dn,,dnz 



ctt + ^cfl + V2e(l + e)o-LT cos 0* 



+e(TTT cos2(?!)* +/ie^2e(l -e)(T£T' sin0* . (5) 

Here = [a/ {2t:''Q'')]{E,, / E,)[\q:^\/ [l - e)]- e is de- 
fined by the electron scattering angle 9^ and the ^ 

2 Analysis and Results [ 

We parameterize the bare "f* N N* vertex func- 
tions r^. . ;,„((7,Q') as 



EnIq) V ko| 
xG,(7V*,Q2)(5,,(,^_,„), (6) 

where qu and qo are defined by Mjv = qR + EpfiqR) 
with N* mass and W = qo+Ej^{qo), respectively, and 

GxiN'.^Q^) = AxiN^Q'^), transverse photon(7) 
= SxiN*,Q'^), longitudinal phot^a) 

For later discussions, we also cast the helicity ampli- 
tudes of the dressed vertex Eq. dU into the form of 



and p{e,e'TT^ 



up to M/ = 1.6 GeV and 



1.45 (GeV/c)^. The availability of the data in the 
corresponding {W,Q^) region are found in Table [TJ 
The resulting parameters are then confirmed against 
the original five- fold differential cross section data ^ . 
This procedure could overestimate/underestimate the 
errors of our analysis, but is sufficient for the present 
exploratory investigation. 

In section II, we present the results from our anal- 
ysis. Discussions on future developments are given in 
section HI. 



Table 1. values, in GeV, for which there are 
available structure function data (Q^ < 1.45 
(GeV/c)2). Data from Refs. 0, [l3,lli, [H 







0"LT' 






0.4, 0.525, 0.65, 0.75, 0.9, 1.15, 1.46 


0.4,0.65 


7*p- 




0.3, 0.4, 0.5, 0.6 


0.4,0.65 



Eq. ([6]) with dressed helicity amplitudes 

Ax{N*,Q') = Ax{N',Q') + Ar-{N\Q'), (9) 
Sx{N*,Q') = Sx{N%Q') + Sr-{N%Q'),{lO) 

where A'^-^-^N* ,Q'-^) and S^-^iN* ,Q^) are due to the 
meson cloud effect defined by the second term of 
Eq. Q. 

With the hadronic parameters fixed in analyz- 
ing the TriV reaction data the only freedom in 
analyzing the electromagnetic meson production re- 
actions are the electromagnetic coupling parameters 
of the model. If the parameters listed in Ref. 
are used to calculate the non-resonant interaction 



vi's'M'B',x-,x,^ik',q) in Eqs. Q and the only pa- 
rameters to be determined from the data of pion elec- 
troproduction reactions are the bare helicity ampli- 
tudes defined by Eq. ®. Such a highly constrained 
analysis was performed in Ref. 0] for pion photo- 
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W = 1.14(GeV) W=1.22(GeV) W = 1.30(GeV) W = 1.42 (GeV) W=1.50(GeV) W = 1.58(GeV) 
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Fig. 1. Fit to p{e,e'-!Y°)p structure functions at = 0.4 (GeV/c)^. Here 6 = 9'^. The solid curves are the 

results of Fitl, the dashed curves are of Fit2, and the dotted curves are of Fit3. (See text for the description 

of each fit.) The data are taken from Refs. Ell- 



production. It was found that the available data of 
jp — > 7r°p, TT+n can be fitted reasonably well up to 
invariant mass W <1.6 GeV. In this work we extend 
this effort to analyze the pion electroproduction data 
in the same W region. 

We first try to fix the bare hclicity amplitudes 
by fitting to the data of <7t + l , o'lt , and Utt of 
p{e,e'TT°)p in Ref. which covers almost all {W,Q^) 
region we are considering (see Table. [TJ . In a purely 
phenomenological approach, we first vary all of the 
helicity amplitudes of 16 bare TV* states, considered 
in analyzing the nN ttN, ttttN data , in the fits 
to the data. It turns out that only the helicity ampli- 
tudes of the first A''* states in 5*11 , Pn , P33 and D^^ are 
relevant in the considered W < 1.6 GeV. Thus only 
the bare helicity amplitudes associated with those 
four bare TV* states (total 10 parameters) are varied 
in the fit and the other bare helicity amplitudes are 
set to zero. The numerical fit is performed at each 

independently, using the MINUIT library. 

The results of our fits are the solid curves in the 
top three rows of Fig.[U Clearly our results from this 
fit agree with the data well. We obtain similar qual- 
ity of fits to the data of Ref. 10| at other values 
listed in Table. [T] We have also used the magnetic 
All form factor of 7*A^ A(1232) extracted from 
previous analyses as data for fitting. We refer the 
results of this fit to as "Fitl" . 



The helicity amplitudes of Su, Pn, and D13 re- 
sulting from Fitl are shown in Fig. [2j The solid cir- 
cles are the absolute magnitude of the dressed he- 
licity amplitudes (O and (fTU|) . The errors there are 
assigned by MIGRAD in the MINUIT library. More 
detailed analysis of the errors is perhaps needed, but 
will not be addressed here. The meson cloud ef- 
fect (dashed curves), as defined by and S'™''- 
of Eqs. (O and (|10p and calculated from the second 
term of Eq. (U) , are the necessary consequence of the 
unitarity conditions. They do not include the bare he- 
licity term determined here and are already fixed in 
the photoproduction analysis ^. Within our model 
(and within Fitl), the meson cloud contribution is 
relatively small in 5*11 and ^1/2 of Dyi, even in the 
low region. 

Here we note that our helicity amplitudes defined 
in Eqs. (O and pUjl are different from the commonly 
used convention, say v4™" and 5*™^, which are ob- 
tained from the imaginary part of the 7*iV — > TriV 
multipole amplitudes '^''l . This definition leads to he- 
licity amplitudes which are real, while our dressed 
amplitudes are complex. It was shown in Ref. [lit 
that for the A(1232) resonance our dressed helicity 
amplitudes ^ and pO|l can be reduced to ^4™^ and 
S"™", if we replace the Green function G^at with its 
principal value in all loop integrals appearing in the 
calculation. However, such reduction is not so triv- 
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Fig. 2. Extracted helicity amplitudes for Sn at = 1535 MeV (upper panels), Pn at W = 1440 MeV (middle 
panels), and D13 at = 1520 MeV (lower panels). Solid points are from Fitl; dashed curves are the meson 
cloud contribution. Open circles and triangles at — OA (GeY/c)^ are from Fit2 and Fit3, respectively. 



ial for higher resonance states because the unstable 
7rA,p7V, (T-/V channels open, and thus the direct com- 
parison of the helicity amplitudes from other analyses 
becomes unclear. 

At ^ 0.4 (GeV/c)2, the data of all struc- 
ture functions both for p(e,e'7r'')p and p(e,e'7r+)n are 
available as seen in Table. [TJ To see the sensitivity 
of the resulting helicity amplitudes to the amount of 
the data included in the fits, we further carry out 
two fits at this Q^, referred to as Fit2 and Fit3, re- 
spectively. Fit2 (Fit3) further includes the data of 
Refs. 



12 



13(](Ref. [lH) in the fit in addition 
to those of Ref. [10| which are used in Fitl. This 
means that Fit2 includes all available data both from 
p{e,e'TT°)p and p(e,e'7r+)n, whereas Fit3 includes the 
same data but from p{e,e'TT°)p only. The results of 

It has been seen in Fig. [3] that all of our current 
fits underestimate CTt of p(e,e'7r+)ri at forward angles. 
We find that this can be improved by further vary- 
ing the S31 and P13 bare helicity amplitudes within a 
reasonable range. We confirm that the same conse- 
quence is obtained also at other W, and find that the 
P13 has contributions mainly at low (high) W. 

We also find that the inclusion of the bare S31 and 
P13 helicity amplitudes does not change other struc- 
ture functions than ctt of p(e, e'7r"'")n (the change is 
within the error). This indicates that those two he- 
licity amplitudes are rather relevant to p(e,eV+)n, 
but not to p{e,e'n'^)p. As shown in Table (TJ however, 
no enough data is currently available for p{e,e''K'^)n 
above = 0.4 {GeV/cf. The data both of the 
p{e,e'n°)p and p(e,e'7r+)n at same values are de- 
sirable to pin down the dependence of the 6*31 and 
P13 helicity amplitudes. 



each fit are the dashed and dotted curves in Fig.[l]for 
p{e,e'TT°)p and Fig. [3] for p(e,e'7r+)n, respectively. 

The corresponding change in the dressed helicity 
amplitudes are also shown as open circles and trian- 
gles in Fig. [21 A significant change among the three 
different fits is observed in most of the results. This 
indicates that fitting the data listed in Table [T] are far 
from sufficient to pin down the 'y* N —^ N* transition 
form factors up to — 1.45 (GeV/c)^. It clearly indi- 
cates the importance of obtaining data from complete 
or over-complete measurements of most, if not all, of 
the independent p{e,e''K)N polarization observables. 
Such measurements were made by Kelly et al. ^ in 
the A (1232) region and will be performed at JLab 
for wide ranges of W and in the next few years ^ . 



We now turn to show the coupled-channels ef- 
fects (CCE). In Fig. m we see that when only the 
ttN intermediate state is kept in the M'B' summa- 
tion of the non-resonant amplitude [Eq. ^] and the 
dressed 7*A^ — > N* vertices [Eq. ([U], the predicted 
total transverse and longitudinal cross sections <Tt 
and <7l of p{e,e'Tr°)p are changed from the solid to 
dashed curves. This corresponds to only examining 
the CCE on the electromagnetic (Q^-dependent) part 
in the 7*A^ ^ ttN amplitude. All CCE on the non- 
electromagnetic interactions are kept in the calcula- 
tions. We find that the CCE tends to decrease when 

increases. This is rather clearly seen in CTt- In 
particular, the CCE on <Tt at high W ^ 1.5 GeV is 
smaU (10-20%) aheady at = 0.4 (GeV/c)^. (The 
effect is about 30-40% at = 0.) This is understood 
as follows. In Eq. ^ we can further split the reso- 
nant amplitude as ^barc+^^.c where t^^^^ and 
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-0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 




-0.5 0.5 


-0.5 0.5 


-0.5 0.5 


-0.5 0.5 


-0.5 0.5 


-0.5 0.5 












^ — 






■ ■ • 

, 1 , I'.il, , 


i t\ " :^l - ," 1 , 





I I I.I. I'. ill. I iT V l-V l I 1-7- 

-0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 
cos9 cosO cosO cosO cos6 cosO 



Fig. 3. Structure functions of p(e, e'Tv^)n at — 0.4 (GeV/c)^. Here 9 = 9'!^. The solid curves are the resuhs 
of Fitl, the dashed curves are of Fit2, and the dotted curves are of Fit3. (See text for the description of 
each fit.) As for the a^r' , results at W = 1.14, 1.22, 1.3, 1.38, 1.5, 1.58 GeV (from left to right of the bottom 
row) are shown, in which the data are available. The data in the figure are taken from Ref. [l^. [l3|. 



*2.c. 9.re the same as but replacing fj^r.^A^Ajv ^^^^ 
its bare part Tj^, ,^ and meson cloud part [the sec- 
ond term of Eq. respectively. The CCE shown 
in Fig.Hcomes from tis^^M,x^Xj^ and t^.,.. We have 
found that the relative importance of the CCE in each 
part remains the same for increasing Q^. However, 
the contribution of non-resonant mechanisms both 
on ^isjvTTiv.A Ajv and t^-^ ^ to the structure functions 
decreases for higher compared with itarc- This 
explains the smaller CCE compared with the photo- 
production reactions ^ . The decreasing non-resonant 
interaction at higher is due to its long range na- 
ture, thus indicating that higher reactions provide 
a clearer probe of N* . We obtain similar results also 
for p(e,e'7r+)n. 

It is noted, however, that the above argument does 
not mean CCE is negligible in the full 7*iV irN re- 
action process. In the above analysis we kept the 
CCE on the hadronic non-resonant amplitudes, the 
strong N* vertices, and the N* self-energy, which are 
Q^-independent and remain important irrespective of 
Q^. We have found in the previous analyses that 
the CCE on them is significant in all energy region 
up to IV = 2 GeV. 



3 Summary and outlook 

We have explored how the available p(e,e'7r)iV 
data can be used to determine the 7*iV — > TV* 
transition form factors within a dynamical coupled- 
channels model. Within the available data, the 
7*iV N* bare helicity amplitudes of the first N* 
states in Sn, Pn, P33 and Di^ can be determined in 
the considered energy region, W < 1.6 GeV. We fur- 
ther observe that some of these parameters can not 
be determined well. The uncertainties could be due 
to the limitation that only data of 4 out of 11 inde- 
pendent p{e,e''K)N observables are available for our 
analysis. The data from the forthcoming measure- 
ments of double and triple polarization observables 
at JLab will be highly desirable to make progress. 

We found that the underestimation of the ar 
of p{e,e' 1:^)71 at forward angles can be improved 
by further considering the S-ix and P13 bare helic- 
ity amplitudes. Furthermore, these amplitudes can 
have relevant contribution to p{e,e'Tr'^)n, but not to 
p(e,e'7r°)p. The p(e,e'7r+)n data of wide region as 
well as p(e,e'7r+)n seem necessary for determining the 

dependence of the 5*31 and P13 helicity amplitudes. 

For testing theoretical predictions from hadron 
structure calculations such as LQCD, the quantities 
of interest are the residues of the ^ ttN ampli- 
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Fig. 4. Coupled-channels effect on the integrated structure functions ariW) and (7_l(W) for = 1.45 
(GeV/c)^ for p{e,e'Tv'^)p reactions. The solid curves are the full results calculated with the bare helicity 
amplitudes of Fitl. The dashed curves are the same as solid curves but only the nN loop is taken in the 
M'B' summation in Eqs. ^ and (g]). 



tudes, defined by Eqs. UJ-llll), at the corresponding 
resonance poles. If the resonance poles are associ- 
ated with the amplitude <f/„^iv,A.,A„(fc,9, of 
Eq. ([3]), the extracted residues are directly related to 
the dressed form factors f ■^.^^,g,J^^,g, (fc', ly). An an- 
alytic continuation method for extracting these infor- 
mation has been developed ^ , and our results along 
with other hadronic properties associated to nucleon 
resonances will be published elsewhere. Let us men- 
tion that the extracted form factors are complex num- 
bers and some investigations are needed to see how 
they can be compared to the usual helicity ampli- 
tudes, which are real numbers, listed by PDG 
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